Renormalization of the Coulomb interaction in one-dimensional electron systems 
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Renormalization group methods are used to study the low-energy behavior of the unscreened 
OO ' Coulomb interaction in a one-dimensional electron system. By applying a GW approximation, a 

0^ ' strong wavefunction renormalization is found in the model, as well as renormalization of the Fermi 

velocity. This may have significant effects depending on the filling level of the system. The effective 
' ' couplings remain bounded at arbitrarily low energies so that, in spite of the long-range character of 

, the interaction, the system still falls into the Luttinger liquid universality class. 
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. The recent experimental availability of single fuUerene nanotubes has renewed the interest in the study of one- 
I— — i' dimensional electron systems 0. In one spatial dimension the Luttinger liquid concept replaces the Fermi liquid 
^ : picture, and provides the paradigm of a system with strong electronic corrrelations. While the fuUerene tubules are 
genuine one-dimensional systems and, therefore, appropriate to test the consequences of the electronic correlations, it 
is also true for the same reason that the bare interaction is not screened in the case of isolated tubules. The analysis 
of the Luttinger liquid behavior is usually made under the assumption of a local or short-range interaction |^-|^ , but 
, the effects that the unscreened Coulomb interaction may have on the properties of the one-dimensional system have 
^ ' not been completely assessed. 

The bosonization approach, in particular, is one of the most elegant and rigorous ways of describing the Luttinger 
liquid behavior. However, the Coulomb interaction is singular at small momentum transfer, and the usual formulas 
that give the renormalized parameters like the Fermi velocity vp m terms of the interactions 52(9) and 54(9) cannot be 
^ ' applied at arbitrarily small energy scales, unless the Coulomb potential is screened by external charges The more 
physical reason for the unfeasibility of the bosonization approach to the present problem regards the construction 
of the ground state of the system in terms of the density operators of the noninteracting theory. Such state is 
normalizable provided that at large momenta ^g2{q) / {2ttv p + 54(9)) ^ while it is unclear that this condition 
may be satisfied in the case of the Coulomb interaction. 

In this Letter we are going to use renormalization group (RG) methods to find the low-energy effective theory of 
the l/\x\ interaction. The main point that we want to address is the stability of the Coulomb interaction in the 
RG framework, and whether it falls into the Luttinger liquid universality class. In order to deal conveniently with 
this question one has to introduce the dynamical screening due to plasmons, what we accomplish by implementing a 
^\ , GW approximation. The suitability of this approximation in the study of one-dimensional systems has been recently 
shown in Ref . . The same approach has been also tested in the study of the crossover from Fermi liquid to Luttinger 
liquid behavior , as well as in the study of singular interactions in dimension 1 < d < 2 . 

To be more precise, we want to pose the problem of a one-dimensional one-band model with an interaction hamil- 
tonian 
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O . H,nt = ^ dxdx' ^+{x)^{x)- -^>+{x')^{x') (1) 
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where ^'(x) is the electron annihilation field operator. The RG method is a sensible approach to deal with this problem 
. since the '^l\x\ interaction potential (as well as the 5(x) potential) gives rise to a marginal four-fermion interaction. 
?H The scaling dimension of the electron field "^{x) is —1/2 , in length units. This means that the interaction hamiltonian 
. 5^ . in (|l|) scales appropriately, with a dimensionless coupling constant (in units in which ?i = c = 1), as the energy 
scale is reduced down to the Fermi level flOl . 

However, the drawback in dealing with (|l|) is that it contains a highly nonlocal operator, which makes unclear the 
applicability of RG methods, usually devised to deal with a set of local operators. This problem can be circumvented 
by introducing a local auxiliary field to propagate the Coulomb interaction. The interaction potential in real space 
l/|a; — x'l corresponds to a singular interaction ~ log \q\ in momentum space. Therefore, it makes sense to focus on 
the electron states that lie close to any of the two Fermi points of the band at -izkp , and to consider two respective 
left (L) and right (R) branches with linear dispersion relation for such relevant states. In practice, one enforces 
this approximation by considering the excitations that have less energy than a given bandwidth cutoff Ec . The 
hamiltonian can be written in the form 
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H = iVF J dx {^+{x)d^^R{x) - ^+ix)d.,^L{x)) +ej dx {^+{x)-^l{x) + (bix) (2) 

where the (j){x) field propagates the interaction 

i{mx, t) t')) = ^5{t - t')^^ (3) 

In writing (H) we have neglected backscattering processes that connect the two branches of the dispersion relation. 
This is justified, in a first approximation, as for the Coulomb interaction the processes with small momentum transfer 
have much more strength than those with momentum transfer ^ 2kp . The backscattering processes give rise, however, 
to a marginal interaction, that should be taken into account after renormalization of the dominant forward scattering 
channel. 

We may think of (j){x) as the scalar potential in three-dimensional quantum electrodynamics. However, the differ- 
ences with that theory in the present case are notorious since the propagation of (j){x) is that of a genuine field in three 
spatial dimensions, while the electrons are confined to one dimension. In general, one may expect a better infrared 
behavior in the present model. The propagator of the (j}{x) field in momentum space can be read from the relativistic 
expression , after sending the speed of light to infinity. 

The usual one-dimensional propagator ^ logdgj/A) is recovered from (^) upon integration of the dummy variables 
Qy and Qz- We remark that the ultraviolet cutoff A for excitations along the y and z transverse directions is needed 
when projecting the three-dimensional interaction down to the one-dimensional system. 

The usefulness of the representation (Q) can be appreciated in the renormalization of the model at the one-loop 
level. We study the scaling behavior of the irreducible functions as the bandwidth cutoff Ec is sent towards the Fermi 
level, Ec ^ 0. The self-energy to the one-loop order is 



2t: J^^ 2t: + v^p + kf J (2^)2 + pi + p- 



The limit A; — > has to be taken carefully in this expression, by first combining the two denominators with the use of 
Feynman parameters . Finally we get 
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The term linear in k in l](fc,0) represents a renormalization of the Fermi velocity, which grows upon integration of 
the high-energy modes. Obviously, there is no correction linear in lo^ renormalizing the electron wavefunction at the 
one-loop level. This is consistent with the fact that the integration of high-energy modes at ~ Ec does not renormalize 
the three-point vertex F . We want to stress the difference of the logarithmic renormalization of Vp with respect to 
the usual finite corrections due to a short-range interaction. The nontrivial scaling of vp is a genuine effect of the 
long-range Coulomb interaction, which also takes place in higher dimensions [T^ ]. 

Incidentally, the above computation exemplifies how the Ward identity that ensures the integrability of the Luttinger 
model does not hold in the present case. The Ward identity is a relation between the electron Green function G{p, ujp) 
and the three-point vertex T{p,LUp; fc, Wfc) at a given branch For the right-handed modes, for instance, it is 

J., , ^ G'^fa^p) - G-\p-k,ujp-u;k) 

T{p, ujp] k, uJk) = ; (7) 

UJk — Vpk 

By focusing on the singular dependences on the bandwidth cutoff E^ one can check that is already violated in 
our model to first order in perturbation theory. Actually, the dependence of the vertex F on the variables {k,ujk) of 
the external interaction line is 

■T-^ I. ^ ^ ■ fc log A 

ir(p,ujp;k,ujk) ^ -i-^ r (8j 

47r^ LiJk — Vpk 



2 



We notice that the Ward identity would be satisfied if the scaling could be implemented simultaneously in the 
transverse ultraviolet cutoff A and the bandwidth cutoff E^- However, even in a real system the scaling in A gets 
locked by the finite cross section of the wire, while only the scaling in the longitudinal direction operated by Ec is 
allowed. In this respect, the gauge invariance of quantum electrodynamics is broken by the anisotropy of the electron 
system, as felt by the propagation of the three-dimensional electromagnetic field. 

The renormalization of vp at the one-loop level is not, in general, a sensible effect from the physical point of view, 
since the propagator of the 4'{x) field is drastically modified by the quantum corrections. Its one-loop self-energy is 
given by the sum of two respective particle-hole diagrams with modes in the left and the right branch of the dispersion 
relation. The result for both fermion loops is finite, and they combine to give 

zn(fc,a;,)=z-^,^,_^, (9) 

Taking into account these particle-hole processes leads to a modified propagator of the (l){x) field 

*(</)(fc, c.) ^(-fc, -u;)) = 1/ (- ^^J^y^^ + n(fc, (10) 

The expression ( p^ ) provides a sensible approximation for the scalar propagator, as it takes into account the plasmons 
in the model. Quite notably, the result ( p^ ) turns out to be the exact propagator in the Luttinger model with 



short-range interactions |13|. 

Our approach is that of using the scalar propagator in the renormalization of the Fermi velocity and the 
electron wavefunction. For this purpose we compute the electron self-energy by replacing the interaction in Fig. |l]by 
the dressed interaction (pO|) 



2m{^j,+u:k)-VF{p + k)i.^-^^logi\p\/A) 



This corresponds to taking the leading order in a expansion, in a model with N different electron flavors. In our 
case, such approximation to the self-energy is also justified since it takes into account, at each level in perturbation 
theory, the most singular contribution at small momentum transfer of the interaction. Due to the cancellation of 
fermion loops with more than two interaction vertices which still takes place in the same way as in the Luttinger 
model, the representation ( pT| ) for the self-energy only misses the effects of vertex and self-energy corrrections, but 
these are treated consistently in the RG framework. 



The only contributions in (11) depending on the bandwidth cutoff are terms linear in ujk and k. In this respect, it is 
worth mentioning that, although the usual perturbative approach gives rise to poles of the form k"^ /{ujk — vpk) in the 
self-energy p^, these do not arise in the GW approximation. There is no infrared catastrophe at uJk ~ vpk, because 
of the correction in the slope of the plasmon dispersion relation with respect to its bare value vp- The contribution at 
small k has to be dealt with particular care, anyhow, as the interaction is not invariant under shifts in the momentum 
transfer. The result that we get for the renormalized electron propagator is 

G^^{k,ujk) = Z^'^ {ujk - Vpk) - Y.{k,ujk) 

^Z, i.,-vpk)-Z, i.,~vpk)J 



/-^-rfp y7b)-4/3 + l/(3/(p)3/2) 
47^' J \P\ (1 - f{p)f 

where f{p) = 1 — log(|p|/A)/(27r^t;i?) and Z^f^ represents the scale of the bare electron field compared to that of 
the cutoff-independent electron field 

■^bareiEa) - 4^'* (13) 

In the RG approach, we require the cutoff-independence of the renormalized Green function, since this object leads 
to observable quantities in the quantum theory. For this purpose, the quantities Z^j, and vp have to be promoted to 
cutoff-dependent effective parameters, that reflect the behavior of the quantum theory as i?c and more states are 
integrated out from high-energy shells of the band. We get then the RG flow equations 
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One can check that Eq. (|l5|) leads to an enhancement of the Fermi velocity at low energies, which goes in the 
direction of screening the effective interaction. Eq. ( p^ ) also corresponds to a sensible effect, as the leading behavior 
is that of suppressing the electron quasiparticle weight. It can be shown that, within the same level of the 1/iV 
approximation, the three-point vertex only gets the cutoff dependence given by the wavefunction renormalization in 
(|l|). This means that the electron charge is not renormalized at low energies in our local field theory framework. 
The behavior of the effective interaction is, therefore, completely encoded in Eq. (p^. 

The RG equation for the effective coupling constant g = e"^ / (Att'^vf) is 

- 4(k^ - I + ^k^) ^''^ 

This equation actually controls all the physical properties of the electron system and, in particular, the vanishing 
of the electron quasiparticle weight. In our model the effective coupling constant g displays marginal behavior. 
However, the logarithmic corrections to scaling in ( ]l6| ) reduce the flow in the infrared. That equation could still 
admit a solution of the form g ~ —go/ log Ec, but this is not realized in the present model as the equation 2gQ = 
+ (70 — 4/3 + 1/(3 1/(1 + 50)"^) does not have any real solution. The flow of g{Ec) given by ( p^ ) seems to be arrested 
close to some fixed-point value, after which it becomes insensitive to further scaling in the infrared. A plot of the flow 
for different values of the bare coupling constant is given in Fig. ^ . 

We see therefore that the Coulomb interaction remains long-ranged in the low-energy effective theory, while the 
effective coupling has a stable flow in the infrared. This means that the one-dimensional system falls into the Luttinger 
liquid universality class. It can be seen from ([l^, in particular, that the electron quasiparticle weight vanishes in the 
low-energy effective theory. The fast wavefunction renormalization shown in Fig. ^ seems to imply a non-algebraic 
behavior of the electron propagator with regard to the frequency and momentum dependence, what is reminiscent of 
similar features found for some correlation functions in the bosonization approach |l4| . 

Moving to the possible implications for real systems, the degree of renormalization of vp may depend, in general, on 
the details of the band structure, as the result (|l^) has been obtained by using the straight form of the polarizability 
for all values of the momentum. The issue of the renormalization of wf is important since it implies a reduction 
in the strength of additional interactions in the system, whether short or long-range, as the effective couplings are all 
given by the couplings in the interaction hamiltonian divided by vp- 

In particular, one may envisage real systems where the Fermi level is close to the top or the bottom of the band, so 
that the polarizability significantly deviates from (^. For very low filling, for instance, it vanishes at large momentum 
transfer like kp/Ec, meaning that the RG fiow should be computed according to the renormalization given by (|^), 
until the dynamical screening due to plasmons takes place. We have represented in Fig. ^ the flow of the effective 
coupling under these circumstances. In this situation there may be a significant renormalization ofvp that affects every 
effective interaction in the model. They flow towards a nonvanishing fixed-point, in any event, since at sufficiently low 
energies there has to be a point in which one recovers the Luttinger liquid regime, where the RG flow equations jisl ) 
and ( p^ apply. The signiflcant renormalization of may also be present in small chains, contributing to explain 
the insulator-metal transition by the effect of the Coulomb interaction observed in the exact diagonalization of finite 
rings 

In conclusion, we have treated the long-range Coulomb interaction in a one-dimensional electron system by applying 
RG methods on a 1 /N truncation of the theory. We have seen that the Coulomb interaction remains long-ranged in 
the low-energy effective theory, and that this feature is compatible with Luttinger liquid behavior of the system. This 
does not seem to support any phase transition with regard to the strength of the interaction. However, significant 
renormalizations of all the effective interactions may arise depending on the band structure, what could be a relevant 
matter in developping the phenomenology of real one-dimensional electron systems and, in particular, the fuUerene 
nanotubes [pi. 
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FIG. 1. Diagram contributing to the electron self-energy. The double line stands for the interaction propagator computed 
in the RPA, according to Eq. (10) in the text. 



FIG. 2. Flow of the effective coupling constant for different bare values, computed in GW approximation (full lines) and 
without the dynamical screening of plasmons (dashed lines). 



FIG. 3. Wavefunction renormalization for bare coupling constant g = 5.0 (thick line) and g = 1.0 (thin line). 
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